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Using a sample of 58 million J/t/j events collected with the BESII detector at the BEPC, more 
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than 100,000 J/^jj — > ppn^ events are selected, and a detailed partial wave analysis is performed. 
The branching fraction is determined to be Br{J/ip ppTr") = (1.33 ± 0.02 ± 0.11) x 10"=^. A 
long-sought 'missing' A'^*, first observed in J/tp pnn~ , is observed in this decay too, with mass 
and width of 2040^4 ± 25 MeV/c^ and 230t| ± 52 MeV/c^ respectively. Its spin-parity favors f"*". 
The masses, widths, and spin-parities of other N* states are obtained as well. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.20.Gk, 14.40.Cs 



I. INTRODUCTION 

Studies of mesons and searches for glueballs, hybrids, 
and multiquark states have been active fields of research 
since the early days of elementary particle physics. How- 
ever, our knowledge of baryon spectroscopy has been 
poor due to the complexity of the three quark system 
and the large number of states expected. 

As pointed out by N. Isgur [J] in 2000, nucleons are 
the basic building blocks of our world and the simplest 
system in which the three colors of QCD neutralize into 
colorless objects and the essential non-abelian charac- 
ter of QCD is manifest, while baryons are sufficiently 
complex to reveal physics hidden from us in the mesons. 
The understanding of the internal quark-gluon structure 
of baryons is one of the most important tasks in both 
particle and nuclear physics, and the systematic study 
of baryon spectroscopy, including production and decay 
rates, will provide important information in understand- 
ing the nature of QCD in the confinement domain. 

In recent years, interest in baryon spectroscopy has 
revived. For heavy baryons containing a charm or bot- 
tom quark, new exciting results have been obtained since 
the experimental evidence for the first charmed baryon 
E+"'" was reported by BNL [^l in 1975 in the reac- 
tion Vfj^p — > //^ATT+TT^Tr+TT^. Many charmed baryons 
have been observed in recent years in CLEO, the two 
B-factories, the Fermilab photo-production experiment, 
FOCUS, and SELEX I, i, MMM- Only a few baryons 
with beauty have been discovered so far. Earlier results 
on beauty baryons were from CERN ISR and LEP [8] 
experiments, while new beauty baryons are from CDF 
and DO at the Tevatron 0, M, Most information 

on light-quark baryons comes from ttN or KN elastic 
or charge exchange scattering, but new results are being 
added from photo- and electro-production experiments 
at JLab and the ELSA, GRAAL, SPRINGS, and MAMI 
experiments, as well as J/ip and ipi'^S) decays at BES. 
However, up to now, the available experimental informa- 
tion is still inadequate and our knowledge on N* res- 
onances is poor. Even for the well-established lowest 
excited states, -/V(1440), N{1535), etc., their properties, 
such as masses, widths, decay branching fractions, and 



spin-parity assignments, still have large experimental un- 
certainties [12]. Another outstanding problem is that, the 
quark model predicts a substantial number of N* states 
around 2.0 GeV/c^ [H, El but some of these, the 
'missing' N* states, have not been observed experimen- 
tally. 

J/ip decays provide a good laboratory for studying not 
only excited baryon states, but also excited hyperons, 
such as A*, S*, and S* states. All N* decay channels 
which are presently under investigation in photo- and 
electro-production experiments can also be studied in 
J/ip decays. Furthermore, for J/ip ^ NNtt and NNtttt 
decays, the Ntt{Ntt) and NTnT{NTnr) systems are ex- 
pected to be dominantly isospin 1/2 due to that the 
isospin conserving three-gluon annihilation of the con- 
stituent c-quarks dominates over the isospin violating de- 
cays via intermediate photon for the baronic final states. 
This makes the study of N* resonances from J/4> de- 
cays less complicated, compared with ttA^ and exper- 
iments which have states that are a mixture of isospin 
1/2 and 3/2. 

N* production in J/ip —^ pprj was studied using a 
partial wave analysis (PWA) with 7.8 x l^fj/ip BESI 
events [16] . Two N* resonances were observed with 
masses and widths of M = 1530 ± 10 MeV, T = 95 ± 25 
MeV and M = 1647 ± 20MeV, T = 1451^^ MeV, and 

spin-parities favoring J^_= ^ • In a recent analysis 
of J/tp — > pnTT^ + c.c. [17], a 'missing' N* at around 
2.0 GeV/c^ named iV2:(2065) was observed, based on 
5.8 X lO"^, J/tp events collected with BESII at the Beijing 
Electron Positron Collider (BEPC). The mass and width 
for this state are determined to be 2065 ± 3^30 MeV/c^ 
and 175 ± 12 ± 40 MeV/c^, respectively, from a simple 
Breit-Wigner fit. In this paper, the results of a partial 
wave analysis of J/ip ppn^ are presented, based on the 
same event sample. 



II. DETECTOR AND DATA SAMPLES 

The upgraded Beijing Spectrometer detector, is a large 
solid-angle magnetic spectrometer which is described in 
detail in Ref. [l^]. The momenta of charged parti- 
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cles are determined by a 40-layer cylindrical main drift 
chamber(MDC) which has a momentum resolution of 

Ip = 1.78%Vl' + p2 (p in GeV/c). Particle identifica- 
tion is accomplished by specific ionization (dE/dx) mea- 
surements in the drift chamber and time-of-flight (TOF) 
information in a barrel-like array of 48 scintillation coun- 
ters. The dE/dx resolution is (JdE/dx = 8.0%; the TOF 
resolution for Bhabha events is ctof = 180 ps. A 12- 
radiation-length barrel shower counter (BSC) comprised 
of gas tubes interleaved with lead sheets is radially out- 
side of the time-of-flight counters. The BSC measures 
the energy and direction of photons with resolutions of 
as/E ~ 21%/y/E {E in GeV), = 7.9 mrad, and 
CTz = 2.3 cm. Outside of the solenoidal coil, which pro- 
vides a 0.4 Tesla magnetic field over the tracking volume, 
is an iron flux return that is instrumented with three dou- 
ble layers of counters that identify muons of momenta 
greater than 0.5 GeV/c. 

In this analysis, a GEANTS-based Monte Carlo (MC) 
program, with detailed consideration of detector perfor- 
mance is used. The consistency between data and MC 
has been carefully checked in many high-purity physics 
channels, and the agreement is reasonable. More details 
on this comparison can be found in Ref. [l^ . 



III. EVENT SELECTION 

The decay J/^p —> ppw^ with tt^ — s- 77 contains two 
charged tracks and two photons. The first level of event 
selection for J/ip —>■ ppn^ candidate events requires two 
charged tracks with total charge zero. Each charged 
track, reconstructed using MDC information, is required 
to be well fitted to a three-dimensional helix, be in the 
polar angle region | cos 6 mdc\ < 0.8, and have the point 
of closest approach of the track to the beam axis to be 
within 1.5 cm radially and within 15 cm from the center 
of the interaction region along the beam line. More than 
two photons per candidate event are allowed because of 
the possibility of fake photons coming from interactions 
of the charged tracks in the detector, from p annihilation, 
or from electronic noise in the shower counter. A neutral 
cluster is considered to be a photon candidate when the 
energy deposited in the BSC is greater than 50 MeV, the 
angle between the nearest charged tracks and the cluster 
is greater than 10°, and the angle between the cluster de- 
velopment direction in the BSC and the photon emission 
direction is less than 23°. Because of the large number 
of fake photons from p annihilation, we further require 
the angle between the p and the nearest neutral cluster 
be greater than 20° . Figures [T] (a) and (b) show the dis- 



tributions of the angles 9-yp and O-yp between the p or 
p and the nearest neutral cluster for J/ip ppw^ MC 
simulation; most of the fake photons from p annihilation 
accumulate at small angles. 

To identify the proton and antiproton, the combined 
TOF and dE/dx information is used. For each charged 
track in an event, the particle identification (PID) 
Xp/_d(*) is determined using: 

/ .\ E^OEYjieasured 'EOF^xpectedV^J 
XTOFW = 7TT 

CJTOFyi) 

/ .\ dE / dXjyicasured dE / dx f^xpectedij'^ 

XdE/dxV') — 77\ 



(^dE/dxii) 

Xpioi^) = XdE/dx{-i) + Xtof(0> 



where i corresponds to the particle hypothesis. A charged 
track is identified as a proton if X%id ^'^^ proton hy- 
pothesis is less than those for the n or K hypotheses. 
For the channel studied, one charged track must be iden- 
tified as a proton and the other as an antiproton. The 
selected events are subjected to a 4-C kinematic fit under 
the J/ip — *■ J3P77 hypothesis. When there are more than 
two photons in a candidate event, all combinations are 
tried, and the combination with the smallest 4-C fit 
is retained. 

In order to reduce contamination from back-to-back 
decays, such as J/tp pp etc., the angle between 
two charged tracks, 02chrg, is required to be less than 
175°. Figures [2] (a) and (b) show the distributions of 
P^^ for simulated J/ip ^ ppn'^ and J/ip jpp events, 
respectively. Selected data events are shown in Fig- 
ure. [2] (a). Here, the variable P^^ is defined as: P^^ = 

4|P„iissp sin^ ^7/2 where Pmiss is the missing momentum 
in the event determined using the two charged particles, 
and 6-y the angle between Pmiss and the higher energy 
photon. By requiring Pf^ > 0.003 GeV^/c'^, background 
from J/f/; — s- jpp is effectively reduced. 

The 77 invariant mass spectrum after the above se- 
lection criteria is shown in Fig. [31 where tt'^ and 77 sig- 
nals can be seen clearly. To select J/t/i — > ppir^ events, 
\M^y - 0.135| < 0.03 GeV/c^ is required. Figures H (a) 
and (b) show the invariant mass spectra of AIp^o and 
Mp^o , respectively, and clear N* peaks are seen at around 
1.5 GeV/c^ and 1.7 GeV/c^. The Dalitz plot of this de- 
cay is shown in Fig. [SI and some N* bands are also evi- 
dent. Both the mass spectra and Dalitz plot exhibit an 
asymmetry for TOp^o and TOp^o, which is mainly caused 
by different detection efficiencies for the proton and an- 
tiproton. The re-normalized Mp^o and Afp^o invariant 
mass spectra after efficiency corrections are shown as the 
solid histogram and crosses, respectively, in Fig. [6l and 
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the agreement is better. 

Other possible J/ij) ^ ppn^ backgrounds are studied 
using MC simulation and data. Decay channels that have 
similar final states as J/tp — > ppn'^ are simulated, and 
J/ip —> ppn'^TT^ is found to be the main background chan- 
nel. Surviving J/tj: pp'k^tt'^ events, passing all require- 
ments described above, are plotted as black dots in Fig.m 
The invariant mass distribution of this background can 
be described approximately by phase space. The tt^ side- 
band, defined by 0.2 < {M^^-Q.Ub) < 0.2278 GeV/c^, is 
used to estimate the background from non-Tr*^ final states, 
such as J/ip IPP, etc.. The circles in Fig. d] show 
the contribution from tt'^ sideband events. In the par- 
tial wave analysis, described below, two kinds of back- 
ground are considered, ir^ sideband background and a 
non-interfering phase space background to account for 
the background from J/ip ppn'^Tr'^ . 



IV. PARTIAL WAVE ANALYSIS 

A partial wave analysis (PWA) is performed to study 
the N* states in this decay. The sequential decay pro- 
cess can be described by J/ip pN*{pN*), N*{N*) 
p7r^(p7r'^). The amplitudes arc constructed using the rel- 
ativistic covariant tensor amplitude formalism [20, 21], 
and the maximum likelihood method is used in the fit. 



A. Introduction to PWA 

The basic procedure for the partial wave analysis is 
the standard maximum likelihood method: 

(1) Construct the amplitude Aj_ for the j-th possible 
partial wave in J/ip — > pNx , Nx —>■ pT^^ or J/'0 — > 
pNx,Nx pTT^ as: 

A, = Kroci^x{BW)xAaecay^x, (1) 

where A-'^rod-x amplitude which describes the pro- 

duction of the intermediate resonance TVx, BWx is the 
Breit-Wigner propagator of Nx^ and Adecay-x is the de- 
cay amplitude of Nx. The corresponding term for the 
Nx is obtained by charge conjugation with a negative 
sign due to negative C-parity of J/i^. 

(2) The total transition probability, lj, for each event 
is obtained from the linear combination of these partial 
wave amplitudes Aj as ijJ — \Y,jCjAj\^ , where the Cj pa- 
rameters are to be determined by fitting the data. 

(3) The differential cross section is given by: 

^ = |I],c,^,P + Ffc„ (2) 

where, is the background function, which includes 
TT*^ sideband background and non-interfering phase space 
background. 

(4) Maximize the following likelihood function InC to ob- 
tain Cj parameters, as well as the masses and widths of 
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FIG. 2: Pt?y distributions. For (a), crosses are data, and the histogram is MC simulation of J/ip 



■ ppn^ . 



simulated J/^ — > jpp events. P^^ = 4|P, 



' sin 6^/2 where Pmis 



(b) Distribution for 
is the missing momentum in the event determined using 



the two charged particles, and 9^ is the angle between Pmiss and the higher energy photon. 
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FIG. 3: The 77 invariant mass spectrum of J 
candidates. 
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the resonances. 



InC = In 

k=l 



w(6 



(3) 



where ^fc is the energy-momentum of the final state of 
the k-ih observed event, w(^) is the probability to gen- 



erate the combination ^, e(f) is the detection efficiency 
for the combination f . As is usually done, rather than 
maximizing C, S = — ln>C is minimized. 

For the construction of partial wave arnplitudes, we 
assume the effective Lagrangian app roach I23I. [23 | with 
the Rarita-Schwinger formalism [20I [2^. [251. l26j . In this 
approach, there are three basic elements for construct- 
ing amplitudes: the spin wave functions for particles, the 
propagators, and the effective vertex couplings. The am- 
plitude can then be written out by Feynman rules for tree 
diagrams. 

For example, for J/i/; -> iViV*(|^) 
7V(ki, si)Af(K2, S2)7r(K3), the amplitude can be con- 
structed as: 



A:- 



+C3Kii/Kia)75u(ki, si)?/-^, (4) 



where m(k2, S2) and v{ki, si) are i-spinor wave functions 
for N and iV, respectively; ip'^ is the spin-1 wave function, 
i.e., the polarization vector for J/ip. The ci, C2, and 
C3 terms correspond to three possible couplings for the 
J/^ NN*{^~^) vertex. They can be taken as constant 



parameters or as smoothly varying vertex form factors. 

3/24 



The spin |^ propagator P^%^ for N*(^~^) is: 
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FIG. 4: The invariant mass spectra of (a) Mp^o and (b) Mp^o for J/ip ^ pp-K^ candidate events, where the circles with error 
bars are the background events estimated from tt*^ sideband events, and the black dots with error bars are those from simulated 
J/^p — > ppTT^Tv'^ events passing the selection criteria. 
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FIG. 5: Dalitz plot of J/if; 



■ ppTv^ candidates. 



FIG. 6: The re-normalized invariant mass spectra of M^^o 
and Mp^Q after correction for detection efficiency, where the 
histogram is Mp^o and the crosses are Mp^o. 



3/2+ Ml. ~p^ + iMn'Tn' ^ i constructed similarly ^ 



2j,Mpf p^^^^ — p''7'^ The possible intermediate resonances are listed in Ta- 

' (5) bleU Of these states, only a few are (well) established 

states, while A'a;(1885) is one of the 'missing' N* states 
with p = K2 + K3. Other partial wave amplitudes can be predicted by the quark model and not yet experimentally 
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observed. Nx{2065) is also a long-sought 'missing' N* , 
which was observed recently by BES [17]. 

For the lowest lying TV* states, iV(1440), iV(1520), 
and Af(1535), Breit-Wigner's with phase space depen- 
dent widths are used. 



BWx{s) 



mT(s) 



m? + imT[s) ' 



(6) 



where s is the invariant mass-squared. The phase space 
dependent widths can be written as p8| : 
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where Bi (q) (Z = 1, 2) is the standard Blatt-Weisskopf 
barrier factor [20| for the decay with orbital angular 
momentum L and p7rAr(s), P-wAis), and Pr;Ar(s) are the 
phase space factors for ttN, ttA, and rjN final states, 
respectively. 



Pxy{s) 



2qxY{s) 



(10) 



qxY 



Vis - {My + MxVKs - {My - MxW) , . 



where X is vr or 77, F is or A, and qxvis) is the 
momentum of X or y in the center-of-mass (CMS) sys- 
tem of XY. For other resonances, constant width Breit- 
Wigner's are used. 

As described in Ref. [2^, the form factors are intro- 
duced to take into account the nuclear structure^ We 
have tried different form factors, given in Ref. :29|, in 



the analysis and find that for J 
factor preferred in fitting is 



resonances, the form 



FxisTT 



N) 



+ (s, 



N 



(12) 



TABLE I: Resonances considered in the PWA analysis. 



Resonance 


Mass(MeV) 


Width(MeV) 




C.L. 


iV(940) 


940 





1 + 

2 


ofT-shell 


iV(1440) 


1440 


350 


1 + 
2 




iV(1520) 


1520 


125 


3- 
2 




iV(1535) 


1535 


150 


1 - 
2 




Af(1650) 


1650 


150 


1 - 
2 




Af(1675) 


1675 


145 


5- 
9 




Af(1680) 


1680 


130 


5 + 
2 




Af(1700) 


1700 


100 


3- 

!+ 


* + * 


Af(1710) 


1710 


100 


1+ 


* + * 


Af(1720) 


1720 


150 


2 




A^^ISSS) 


1885 


160 


3 - 
1+ 


'missing' A^* 


Af(1900) 


1900 


498 


1+ 


** 


Af(2000) 


2000 


300 


1+ 


** 


Af^(2065) 


2065 


150 


2 


'missing' N* 


Af(2080) 


2080 


270 


3- 
2 


** 


Af(2090) 


2090 


300 


1 - 
!+ 


* 


N{2W0) 


2100 


260 


2 


* 



**** Existence is certain, and properties are at least fairly well 
explored. 

*** Existence ranges from very likely to certain, but further 
confirmation is desirable and/or quantum numbers, branch- 
ing fractions, etc. are not well determined. 
** Evidence of existence is only fair. 
* Evidence of existence is poor. 



where s 



or 



N is the invariant mass squared of TV, tt, and for 
I states, the preferred form factor is 



FN{SnN) = e 



(13) 



'AT* 



Therefore, the above form factors are used in this analy- 
sis. 

In the log likelihood calculation, 7r° sideband back- 
ground events are given negative weights; the sideband 
events then cancel background in the selected candidate 
sample. The J/^p ppn^TT^ background is described by 
a non-interfering phase space term, and the amount of 
this background is floated in the fit. 



B. PWA results 

Well established states, such as Ar(1440), iV(1520), 
A^(1535), iV(1650), 7V(1675), 7V(1680) are included in 
this partial wave analysis. According to the framework 
of soft TT meson theory [g^], the off-shell decay process is 
also needed in this decay, and therefore iV(940) {M — 940 
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P P 



FIG. 7: Feynman diagrams of J/tp pp-K^ for the off-shell 
decay process. 

MeV/c^, r = 0.0 MeV/c2) is also included. Fig. [7| shows 
the Feynman diagram for this process. 



1. Resonances in the 1.7 GeV/c? mass region 

In the M = 1.7 GeV/c^ mass region, three resonances 
iV(1700)(f"), iV(1710)(i'^), and iV(1720)(|"^) [H are 
supposed to decay into pn^pTr) final states. According to 
the Particle Data Group (PDG08) [H, only iV(1720) is 
a well established state. We now study whether these 
three states are needed in J/ip ppn^ . This is in- 
vestigated for two cases, first assuming no N* states in 
the high mass region (> 1.8 GeV/c^), and second as- 
suming iV^(2065), 7V(2080), and Ar(2100) states in the 
high mass region. With no N* states in the M > 1.8 
GeV/c^ mass region, the PWA shows that the signifi- 
cances of iV(1700) and 7V(1720) are 3.2a (AS = 11) and 
0.8cr (AS' = 3), and their fractions are 0.3% and 6%, re- 
spectively; only iV(1710) is significant. When Nr,{206^), 
iV(2080), and iV(2100) are included, the N{niQ) makes 
the log likelihood value S better by 65, which corresponds 
to a significance much larger than 5cr. However, neither 
the A^(1700) nor the A^(1720) is significant. We conclude 
that the A^(1710) should be included in the PWA. 



mass region, the log likelihood value improves by 288, 
which corresponds to a significance of greater than ha, 
when a |"^Af^(2065) is added. Thus, the 7V^(2065) is 
definitely needed in this case, and its mass and width are 
optimized to be M = 205716 MeV/c^ and T = 22Qt\\ 
MeV/c2. 

The significance and spin-parity of Afa;(2065) is fur- 
ther checked under the following four hypotheses (A, B, 
C and D) for the high mass resonances. Case A has 
7V(2080) and 7V(2100) included, case B iV(2080) and 
7V(2000), case C 7V(2000), 7V(2080), and iV(2100), and 
case D 7V(2080), 7V(2090), and 7V(2100). The changes of 
the log likelihood values (AS) , the corresponding signifi- 
cances, and the fractions of Nx{20Qh) are listed in Table 

im when a |^ A'a;(2065) is added in the four cases. The 
log likelihood values become better by 58 to 126 when 
TV:, (2065) is included. Therefore, 7V:,(2065) is needed all 
cases. The differences of log likelihood values for differ- 
ent 7Va;(2065) assignments for the four combinations 
are listed in Table IIIII The assignment of = |^ gives 
the best log likelihood value except for the cases where 
there is large interference. Spin-parity of ^ is favored 
for Af^(2065). 



TABLE II: Changes of log likelihood values (AS), the corre- 
sponding significances, and the fractions of A'^i,(2065), when 
A^a;(2065) is added in the four cases. 



Case 


AS 


significance 


fraction (%) 


A 


126 


> 5cr 


23 


B 


158 


> 5cr 


24 


C 


79 


> 5cr 


16 


D 


58 


> 5cr 


22 



2. iV:,(2065) 

The iV2;(2065), a long-sought 'missing' N* predicted 
by the quark model, was observed in J /ij) — s- pn'K~ -f c.c. 
[iTi] with a mass of 2065 ± 3^30 MeV/c^ and a width of 
175 ± 12 ± 40 MeV/c^, determined fro m a simple Breit- 
Wigner fit. We investigate the need for the A'a;(2065) in 
J/^/i ppiT°. Including the iV(1440), 7V(1520), 7V(1535), 
Ar(1650), A^(1675), iV(1680), A^(171G) and the off-sheU 
decay in the PWA fit, different A'j,(2065) spin-parities 
(J^) and different combinations of high mass resonances 
are tried. If there are no other resonances in the high 



TABLE III: Comparison of log likelihood values for different 
assignments for 7Va;(2065). 





1 + 


1 - 


3 + 


3 - 


5 + 

2 — — 


5 - 


A 


85.8 


49.3 


0.0 


-32.2" 


-36.9" 


34.1 


B 


5.0 


68.5 


0.0 


54.3 


-12.1" 


6.3 


C 


98.1 


39.8 


0.0 


85.6 


76.1 


14.4 


D 


44.2 


45.2 


0.0 


25.0 


36.2 


38.0 



" 780% interference between Nx(2mb) and Af(2080) 
" 529% interference between Af(1680) and Af(2000). 
" 860% interference between 7V(1680) and Af(2000). 
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3. Other resonances in high mass region 



In addition to the observed resonances, -/V(2000), 
7V(2080), 7V(2090) and iV(2100), as well as the iV^(2065), 
there is another possible 'missing' N* state, 7Va;(1885), 
which is predicted by theory but not yet observed. 

a) 7V:r(1885) 

In the p{P)tt'^ invariant mass spectrum, shown in 
Fig. m no obvious peak is seen near 1.89 GeV/c^. We 
study whether this state is needed in the partial wave 
analysis for the four cases. The significances are l.Za 
(AS* = 3.0), 3.2(7 (A5 = 8.8), iAa {AS = 9.7), and 
greater than 5a {AS = 28.0) in cases A, B, C, and D, 
respectively, when a A'a;(1885) is included. Thus, the 
statistical significance is larger than 5(j only in case D. 
In our final fit, 7V2;(1885) is not included. However, the 
difference of including and not including it will be taken 
as a systematic error. 

b) Af(2000), iV(2080), Af(2090), and iV(2100) 

We next study whether 7V(2000), iV(2080), 7V(2090) 
and iV(2100) are all significant in the decay. First, we add 
iV(2000), iV(2080), iV(2090), and A^(2100) one at a time 
with iV(940), iV(1440), 7V(1520), A^(1535), A^(1650), 
iV(1675), A^(1680), 7V(1710), and 7V^(2065) already in- 
cluded. The log hkelihood values get better by 28, 137, 
69, and 73, respectively, which indicates the A^(2080) is 
the most significant, all the significances are larger than 
5a. Second, we include iV^(2065) and 7V(2080) in the 
high mass region and add the other three states A^(2000), 
iV(2090), and A^(2100) one at a time. The significances 
of the iV(2100) {AS = 38) and iV(2090) (AS' = 30) are 
much larger than 5a, while A^(2000) is 3.9ct {AS = 14). 
Third, we include N^{2065), iV(2080), and A^(2100) in 
the high mass region and test whether iV(2000) and 
iV(2090) are needed again. The significances are larger 
than 5a {AS = 23) and 2.7cr (AS" = 7), respectively 
when 7V(2000) and 7V(2090) are included. 

Due to the complexity of the high mass N* states and 
the limitation of our data, we are not able to draw firm 
conclusions on the high mass region. In the final fit, we 
include iV^(2065), A^(2080), and 7V(2100) and take the 
differences of with and without iV(2000) and 7V(2090) as 
systematic errors. 



4- The best results up to now 

We summarize the results we have so far: 

(1) For the three resonances in the M — 1.7 GeV/c^ 
mass region (iV(1700), iV(1710), and A^(1720)), only 
TV (1710) is significant. 

(2) The Nx{2065) is definitely needed in all cases, and 
its spin-parity favors ^ . 

(3) iVa;(1885) is not significant and therefore is not 
included in the final fit. 

(4) For other resonances in the high mass region, 
iV(2080) and A^(2100) are both needed in all cases 
tried, but the other two states iV(2000) and iV(2090) are 
not very significant and so are not included in the final fit. 

Therefore, we consider 7V(940), iV(1440), iV(1520), 
7V(1535), 7V(1650), 7V(1675), 7V(1680), 7V(1710), 
TV:, (2065), Af(2080), and 7V(2100) in the fit. 

Table IIVI lists the optimized masses and widths for 
some A^* resonances; the others are fixed to those from 
PDG08. Here, only statistical errors are indicated. The 
fractions of these states are also listed. 

The Mp^o and Mp^o invariant mass spectra and the 
angular distributions after the optimization are shown in 
Figs. M (a) and (b) and Fig. [9l respectively. In Fig. [8] 
and [HI the crosses are data and the histograms are the 
PWA fit projections. The PWA fit reasonably describes 
the data. 



5. A'^i(1885) significance with optimized N* states 

In the analysis above, the A'a;(1885) was not found to 
be significant. Here its significance is redetermined using 
the optimized masses and widths for the A*'s, and it is 
still only 1.2ct {AS — 2.7). Therefore, we have the same 
conclusion: the A^i,(1885) is not needed. 



6. iV(1900) 

In PDG08 [III, there is an N{1900){^^) state near 

^2,(2065) [l2|. Our previous results show that if there 

Q -t- . . . 

is only one ^ state in this region, the mass and width 

are optimized at M = 2057tg MeV/c^ and F = 220j:i2 
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TABLE IV: Optimized masses and widths, as well as fractions. Errors shown are statistical only. 



Resonance 



Mass(MeV/c2 



145517 
1513ti 
1537lg 
16501^ 
I715I2 



Width(MeV/c2) 



316t| 
127+_l 
1351^ 

230l^ 



J' 



"TT- 
2 

3- 
2 

1- 
2 

1_- 

!+ 

1+ 

2 



Fraction (%) 



Af(1440) 
Ar(1520) 
Ar(1535) 
Ar(1650) 
iV(1710) 

iVa;(2065) 



16.37 
7.96 
7.58 
9.06 
25.33 
23.39 



u 
o 



C 




M(p7t ) 



M(p7t ) 



FIG. 8: The pvr'' and pTr" invariant mass spectra after optimization of masses and widths. Plot (a) is M^^a, and plot (b) is 
Mp^o , where the crosses are data and histograms are fit results. 



MeV/c^, which are consistent with those of 7V2;(2065). 
If iV(1900) is also included in this analysis, i.e. there 
are two | states in this region, we find that the second 

I state also has a statistical significance much larger 
than 5cr (AS' = 49). However, the interference between 
iV(1900) and iV^(2065) is about 80%. This analysis does 
not exclude the possibility that there are two ^ states 
in this region. 



7. Search for additional N* and A* resonances 

Besides the contributions from the well-established N* 
resonances, there could be smaller contributions from 
other N* resonances and even A* resonances from isospin 
violating virtual photon production. 

What might be expected for the isospin violating de- 



cay? For the J/tp pp decay, the isospin violat- 
ing fraction can be estimated using the PDG J/tp lep- 
tonic branching fraction and the proton electromagnetic 
form factor Fpiq"^)^ to be B{J/i; ^ 7* ^ pp) = 
B{J/i; X {Fp{M]f^f = 2.4 x 10"^ The to- 

tal J /'ip pp branching fraction is 2.2 x 10^'' 12]. This 
means, the fraction of decays through a virtual photon 
in the J/ip ^* ^ pp decay mode is close to 1.1%. For 
the non-strange channel, the ratio of photon couplings to 
isospin 1 and isospin is 9:1, so the isospin violating part 
is about 1% for this channel. For the J /ip ppn^ decay, 
one would expect a similar isospin violating fraction. 

If we add an extra state with different possible spin- 
parities ( = i , I , I ) in the large mass (1.65 
GeV/c^ to 1.95 GeV/c^) region with widths from 0.05 
GeV/c^ to 0.20 GeV/c^ and re-optimize, we find that no 
additional A^*'s or A*'s with the statistical significance 
of greater than 5ct are required. 
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FIG. 9: Distributions of (a) the cosine of the sum of the p and p momenta, (b) cosine of the momentum of the p7r° system in 
the pp CMS, (c) cosine of the momentum of the pp system in the pTv^ CMS, and (d) cosine of the momentum of pp in the pTr" 
CMS. The crosses are data and histograms are the fit resuhs. 



8. Search for p{2150) 



V. BRANCHING FRACTION OF J/^P pp-n° 



A resonance with mass 2149 MeV/c^ and ~ 1 is 
listed in PDG08 [H] with the decay p(2150) pp. Here, 
we test whether there is evidence for this decay in our 
sample. The significance of this resonance is less than 3cr 
when we vary the width of this state in the fit from 200 
to 660 MeV/c^. Therefore, our data do not require this 
state. Figure (TU] shows the pp invariant mass spectrum, 
and there is no clear structure near 2149 MeV/c^. 



The branching fraction of J/V' — *■ ppi^^ is obtained by 
fitting the 7r° signal (see Fig. 3) with a ir^ shape obtained 
from J/i/) ppT:^ MC simulation and a polynomial back- 
ground. The numbers of fitted signal and background 
events are 11,166 and 691, respectively. The efficiency 
of J /ip ppTT^ is determined to be 13.77% by MC sim- 
ulation with all intermediate N* states being included. 
Figures 111! (a) and (b) show the p and p momentum dis- 
tributions, where the histograms are MC simulation of 
J/t/j ppir'^ using the J^'s and fractions of N* states 
obtained from our partial wave analysis, and the crosses 
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FIG. 10: The pp invariant mass spectra. Crosses are data and histogram is the PWA fit after the optimization of masses and 
widths. 



are data. There is good agreement between data and MC 
simulation. 

The branching fraction is determined to be: 

Br{J/i; pp7r°) = (1.33 ± 0.02 (stat.)) x 10"^ (14) 

which is higher than that in PDG08 ^ ((1.09 ± 0.09) x 
10-3). 

VI. SYSTEMATIC ERRORS 

The systematic errors for the masses and widths of N* 
states mainly originate from the difference between data 
and MC simulation, the influence of the interference be- 
tween iV(2100) and other states, uncertainty of the back- 
ground, the form-factors, and the influence of high mass 
states, as well as the differences when small components 
are included or not. 

(1) Two different MDC wire resolution simulation 
models are used to estimate the systematic error from 
the data/MC difference. 

(2) In this analysis, the interference between 
A^(2100)(i ) and the low mass regions states such as 

iV(940)(i^) and A^(1440)(i^) can be very large, even 
larger than 50%. We fix the fraction of iV(2100) to be less 
than 10% to reduce the interference and then compare its 
impact on other resonances. The biggest differences for 
the masses, widths, and fractions of the other resonances 
between fixing the fraction of Af(2100) and floating its 
fraction are considered as systematic errors. 



(3) Two kinds of backgrounds are considered in the 
partial wave analysis, 7r° sideband and non-interfering 
phase space. We increase the number of background 
events by 10%, and take the changes of the optimized 
masses and widths as systematic errors. 

(4) Equations (12) and (13) are the form factors used 
in this analysis, where A is 2.0 for N* states with = i 
and A is 1.2 for those with — | and |. Other form 
factors have also been tried, however their log likelihood 
values are much worse than those from the form factors 
used here. We also vary the A values from 2.0 and 1.2 to 
1.5. The biggest differences are taken as the form factor 
systematic errors. 

(5) The effect of using different combinations of states 
in the high mass region on the masses and widths of other 
resonances was investigated above (see Table HH) . and the 
differences also taken as systematic errors. 

Table |V] shows the summary of the systematic errors 
for the masses and widths, and the total systematic errors 
are the sum of each source added in quadrature. 

The systematic errors for the branching fraction 
B{J/ip TT^pp) mainly originate from the data/MC 
discrepancy for the tracking efficiency, photon efficiency, 
particle ID efficiency, fitting region used, the background 
uncertainty, and the uncertainty in the number of J /tp 
events. 

(1) The systematic error from MDC tracking and the 
kinematic fit, 2.18%, is estimated by using different MDC 
wire resolution simulation models. 

(2) The photon detection efficiency has been studied 
using J/i/j ^ /OTT [33|- The efficiency difference between 
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TABLE V: Summary of the systematic errors for masses and widths of A'^* resonances (MeV/c^). 



Systematic error 


iV(1440) 


iV(1520) 


iV(1535) 


iV(1650) 


N{17W) 


iV:,(2065) 




AM 


Ar 


AM 


Ar 


AM 


Ar 


AM 


AF 


AM 


AF 


AM AF 


Data/MC comparison 


3 


14 


2 


13 


2 


11 


4 


12 


1 


12 


1 19 


Interference 


12 


25 


2 


23 


3 


22 


25 


15 


15 


2 


10 20 


Background uncertainty 


18 


51 


11 


23 


6 


28 


2 


8 


4 


10 


5 15 


Different form-factors 


12 


25 


2 


5 


8 


1 


3 


5 


15 


22 


20 14 


Different combinations in high mass region 


35 


21 


7 


12 


5 


10 


5 


23 


20 


35 


10 39 


Total 


43 


67 


13 


37 


12 


39 


26 


31 


29 


44 


25 52 



TABLE VL Systematic error for the branching fraction 
B{J/ip — » Tv^pp) from different sources. 



Sys. error sources 



Wire resolution 

Photon efficiency 

Particle ID 

Mass spectrum fitting 

Number of J/^ events 

Total 



taken as 4% in this decay. 

(4) The TT*^ fitting range is changed from 0.04 - 0.3 

Systematic error(%) GeV/c^ to 0.04 - 0.33 GeV/c^, and the difference , 1.28%, 
is taken to be the systematic error from the fitting range. 
To estimate the uncertainty from the background shape, 
we change the background shape from 3rd order polyno- 
mial to other functions. The biggest change, 1.44%, is 
taken as the systematic error. 



2.18 
4.00 
4.00 
1.93 
4.72 
7.93 



data and MC simulation is about 2% for each photon. 
So 4% is taken as the systematic error for two photons 
in this decay. 

(3) A clean J/^/; — > ppvr+Tr^ sample is used to study 
the error from proton identification. The error from the 
proton PID is about 2%. So the total error from PID is 



(5) The total number of J/tp events determined from 
inclusive 4-prong hadrons is (57.70±2.72) x 10^ 31]. The 
uncertainty is 4.72%. 

Table IVII lists the different sources of systematic er- 
rors for the branching fraction of J/ip — > ppi:^ . The 
total systematic error is the sum of each error added in 
quadrature. 
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TABLE VII: Summary of A'^* states optimized results. 



Resonance 



Mass(MeV/c^) 



1455^7 ± 43 
151314 ± 13 
1537+g ± 12 
1650t^ ± 26 

instj ± 29 

2040li^ ± 25 



width(MeV/c^) 



J' 



Sietg ±67 
127^8 ± 37 
135+8 ± 39 
1451^0 ± 31 
ti ±44 
230l| ± 52 



"TT- 

2 

3- 
2 

1_- 

2 

1_- 

!+ 
1+ 

2 



Fraction(%) 



Branching fraction (xlO 



A(1440) 
Ar(1520) 
A^(1535) 
Ar(1650) 
iV(1710) 
iV^(2065) 



95j 



9.74- 
2.38- 
6.83- 
6.89- 
4.17- 
7.11- 



'25.93 
'10.92 
'15.58 
'27.94 
'30.10 
'24.29 



1.33- 
0.34- 
0.92- 
0.91- 
0.54- 
0.91- 



'3.54 
'1.54 
'2.10 
'3.71 
'3.86 
'3.11 



VII. SUMMARY 

Based on 11,166 J/i/; — > ppTr" candidates from 5.8 x 
10^ BESII J/il> events, a partial wave amplitude analysis 
is performed. A long-sought 'missing' N* , which was 
observed first by BESII in J/tp pTm^ ± c.c, is also 
observed in this decay with mass and width of 2040li^±25 
MeV/c2 and 230^8 ±52 MeV/c^, respectively. The mass 
and width obtained here are consistent with those from 
J/f/; — > pnn^ + c.c. within errors. Its spin-parity favors 
1^. The masses and widths of other N* resonances in the 
low mass region are also obtained and listed in Table IVnl 
where the first errors are statistical and the second are 
systematic. The ranges for the fractions of N* states, 
and thus the branching fractions, are given too. From 
this analysis, we find that the fractions of each N* state 
depend largely on the N* 's used in the high mass region, 
the form factors, and Breit-Wigner parameterizations, as 
well as the background. We also determine the J/ip 
ppTT^ branching fraction to be Br{J/ip ppn^) — (1.33± 



0.02±0.11) X 10"'^, where the efficiency used includes the 
intermediate N* and N* states obtained in our partial 
wave analysis. 
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